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Abstract. The asymptotic behavior of a nonlinear diffusive equation obtained in the framework of the k- 
generalized statistical mechanics is studied. The analysis based on the classical Lie symmetry shows that the 
K-Gaussian function is not a scale invariant solution of the generalized diffusive equation. Notwithstanding, 
several numerical simulations, with different initial conditions, show that the solutions asymptotically ap- 
proach to the K-Gaussian function. Simple argument based on a time-dependent transformation performed 
on the related K-generalized Fokker-Planck equation, supports this conclusion. 

PACS. 05.20.Dd Kinetic theory - 05.20.-y Classical statistical mechanics - 05.90.-|-m Other topics in 
statistical physics 



! 1 Introduction 

In the framework of non-equilibrium thermostatistics, ir- 
. reversible processes can be often described by means of 
' Fokker-Planck equations (FPE's) [T] whose time evolu- 
\ tions are characterized by monotonically non-increasing 
. Lyapunov functionals. In a previous work [5], we derived 
' a non-linear FPE in the picture of a generalized statistical 
, mechanics based on the K-entropy [3,4 , and discussed its 
• relation with the associate Lyapunov functional or Breg- 
man type divergence [5]. Based on the monotonic behav- 
ior of the Lyapunov functional, one can show that any 
' initial localized state, which evolves according to the k- 
, generalized FPE, converges to a stationary solution which 
■ minimizes the Lyapunov functional. In addition, in the 
^ linear drift case, the corresponding stationary solution is 
nothing but the K-Gaussian, which is a generalization of 
Gaussian by replacing the standard exponential with its 
K-generalized version ^ . 

As well known, the diffusive equations play a fundamen- 
tal role in the description of several physical phenom- 
ena. In particular, the linear diffusive equation arising in 
the Boltzmann-Gibbs (BG) statistical mechanics has been 
widely studied jointly with its Gaussian self-similar so- 
lution [7]. Similarly, the g-generalized diffusive equation 
arises, in a natural way, in the framework of the Tsal- 
lis statistical physics. This is equivalent to the porous 
medium equation (PME), which is a transport equation 
of gases or fluids in a porous medium. The main proper- 
ties of PME are now well known [5] and its solutions have 
been widely investigated in literature [QlllOj , in particular 



its self-similar solution is given by the so-called Barenblatt 
solution pn] . 

In this work we study a different nonlinear diffusive equa- 
tion obtained in the picture of the K-generalized statistical 
mechanics. As known, a useful method to obtain the so- 
lutions of a partial differential equation is based on the 
study of its Lie symmetries and the related group invari- 
ant solutions [l2lfT3]. We accomplished such analysis for 
the K-diffusive equation considering only the classical Lie 
symmetries whose generators are functions of the inde- 
pendent and the dependent variables. We will leave the 
study of the generalized Lie symmetries, whose genera- 
tors are functions also of the derivatives of the depending 
variables, to a future work. 

In particular, we show that the only group-invariant so- 
lutions for the K-diffusive equation are the kink-like solu- 
tions (a non-normalizable family of self-similar solutions) 
and the traveling wave solutions (physically irrelevant be- 
cause they are divergent). Next, we explore numerically 
the evolutions of a localized initial state, whose spreading 
is governed by the ^-diffusive equation. It is shown that, 
independently from the initial states, the numerical so- 
lutions approach to a shape which is asymptotically well 
fitted by the K-Gaussian distribution, although this one 
is not a group invariant solution of the equation under 
investigation. 

The paper is organized as follows. In the next section 
2 we briefly review the K-generalized thermostatistics and 
its associated non-linear FPE. In section 3, after present- 
ing an alternative derivation of the K-diffusive equation, 
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we classify its classical Lie symmetries and their related 
group invariant solutions. Section 4 deals with the numer- 
ical analysis and the final section 5 is devoted to summary. 



2 ^-generalized thermostatistics 

We briefly review the generalized thermostatistics based 
on the K-entropy [SUl] given by 



p{v) In^^j p{v) dv 



(1) 



where hi^^j (x) is the K-logarithm deflned us 



2k 



(2) 



The K-entropy S^lp] — S^k,[p\ is an extension of the BG 

entropy by means of a real parameter \k\ £ [0, 1]. 

The inverse function of In^.^^ (a;), namely exp^ y{^)i given 

by 



exp^^j (x) 



+ K^x"^ 



(3) 



is called K-exponential. 

The K-logarithm and ^-exponential functions are the build- 
ing blocks of the statistical mechanics based on the k- 
entropy. In particular, the K-exponential is a monotonic 
increasing and convex function, fulfilling the relationship 
exp^^j(— x) = 1/ expj^^j (a;). Moreover, expj_^j(— oo) = 0, 
exp^^j(O) = 1 and exp^^^{oo) — oo, as the ordinary ex- 
ponential function does. For |x| <C 1, the /c-exponential is 
well approximated by the standard exponential, whereas 
for |x| — > oo, it asymptotically approaches to a power-law: 
exp^^^{x) ^ \2kx\^-^/'^. In the k — > limit, both ln^^j(a;) 
and expj^^j {x) reduce to the standard logarithm and expo- 
nential functions, respectively. Accordingly, the K-entropy 
reduces to the BG entropy. 

Maximizing Sk,[p] under the constraints of the linear 
kinetic energy average and the normalization 



Sp{v) 



-f3 



v'^p{v)dv — j / p{v)dv) —0, (4) 



leads to the K-Gaussian function 



P 



(4(^-4) 



(5) 



Here 7 is the constant fixing the normalization of the dis- 
tribution and depends on the Lagrange multiplier (3 which 
controls the dispersion of the distribution. Finally, the pa- 
rameters a and A are K-dependent constants given by 



1 - K 



A = 



(6) 



respectively, and are related each to the other through the 
relation 

a = expj^j(-l/A) . (7) 

The K-entropy and the corresponding statistical me- 
chanics preserve several properties of the BG theory [3]. In 
particular, the Legendre structure of the thermo-statistics 
theory has been investigated in [14! through the introduc- 
tion of several thermodynamic potentials. 
For instance, it is found that the generalized partition 
function, as a function of /3, given by 



satisfies the Legendre relation 



d 



^ln,,(^.) = -C/, 



(8) 



(9) 



and that the K-entropy and the generalized partition func- 
tion are related through the relationship 



where 



U 



- p{v) dv 



(10) 



(11) 



The function [p] in equation ([5]) is defined by 

oo 

1-n\p\ = / pi'") {p{v)) dv , 



(12) 



where the quantity Uj^j(x) = u^_^y{x), defined as 
, , x'^ + x^'^ 

= — ^ — ) 



(13) 



fulfills the following proprieties u^^y{x) — u^^y{l/x) and 
""{k} ('^) = l/-^- 1^ the K ^ limit, it reduces to the unity: 
= 1- 

Both functions u^^y{x) and Ik[p] play an important role 
in the development of a theory based on the K-entropy. 
Starting from the definition of the partition function , 
we can introduce a K-generalization of the free-energy, ac- 
cording to 

i^Kb] = -^ln,., (^k) , (14) 

In this way, it was shown that the K-free-energy satisfies 
the Legendre transformation structures [14] summarized 
by the following relationships 



f4p]^u[p]-^S4p] 



= u 



(15) 
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2.1 K-generalized Fokker-Planck equation 

In a previous work j2] , we have studied the nonhnear FPE 
associated with the K-entropy. It can be introduced as a 
continuity equation for the density field p{v,t) 



with the current density 



dv \6p{v,t) 



(16) 



(17) 



and 



c4p] = u[p]-ds4p] . (18) 

As shown in [5] the quantity [p] is a K-generalized Lya- 
punov functional, which is monotonically non-increasing 
in time. 

The K-generalized FPE reads 
d d / \ &^ / 

(19) 

where Z? is a constant diffusion coefficient. 
We remark that the current density j(u, t) — j'^"^*(u, t) + 
jdiffj-^^^^ is the sum of a linear drift current j'^'^'^*(w, t) = 
V p, which describe the standard Uhlenbeck-Ornstein pro- 
cess and a nonlinear diffusive current, given by j'^^^{v, t) = 
D9[pU(^j(p)]/9w, which reduces to the standard Pick's 
current j^^'^^{v,t) = Ddp/dv in the k ^ limit. Conse- 
quently, in the same limit, equation (fT9|) reduces to the 
standard linear PPE 



d_ 

dt 



^2 



(20) 



Since dCi^/dt < 0, the functional Ck,[p\ is minimized for 
the stationary solutions of equation (fTOjl according to 



min C^[p] = lim C,^[p] 

t — ^OO 



u[p^']-ds4p'']^f4p'''] 



(21) 



It is easy to verify that the stationary solution of the re- 
generalized FPE is equal to the optimal maximal entropy 
distribution p^^{v), i.e. 



p (v) = lim p{v, t) 

t — ^CxD 



(22) 



= -xp,^^(4(, + l^))^pME(,)^ (23) 
with the position D = 1/(3. 

Without the external potential U the drift therm j'^^'f* 
disappear and equation (|19p reduces to a purely K-diffusive 
equation whose explicit form is given by 

= ,24, 

Equation ([24]) is the main subject of our investigation. 



3 ^-generalized diffusive equation 
3.1 A simple physical derivation 

First, let us describe a physical derivation of the re-diffusive 
equation (I24p . In this respect, we recall that the linear dif- 
fusive equation 

' (25) 
can be obtained from the following physical relations: 



dt 



p{v,t) + -^{y p{v,ty) — (continuity equation) 



dP 
dv 



P oc ln(p(i;,t)) 



(Darcy's law) 
(equation of state) 



(26) 



where v denotes the average velocity of the fluid and P 
denotes the pressure of the fluid, which is related to the 
density of the fluid p(v,t) through a suitable equation of 
state P oc ^(p), with <?(x) = ln(x). 

The continuity equation assures the conservation of the 
zero-th moment 



p{v,t)dv = M 



(27) 



which corresponds to the total mass M of the fluid. 

The PME [5] can be derived from the same physical 
relations by posing, in the equation of state, <l>(x) — x^: 



■g^p{v,t) + -^{v p{v,t)) — (continuity equatic 

(Darcy's law) 
(polytropic fluid) 



dP 

P oc p{v,tY 



(28) 



Here, v is a, positive real parameter characterizing the de- 
gree of porosity of the medium and characterizing also the 
power-law dependence of polytropic fluid. In other words, 

is a polytropic index. 
The corresponding PME is given by 



(29) 



with q = ly + 1. 

Interestingly, the same settings can be also applied to the 
K-diffusive equation (|24p , merely replacing, in the equation 
of state, the logarithm with its ^-generalization: (l>(x) — 
ln{,j(a;). 

In this way, the re-diffusive equation (|24p can be derived 
according to 

■§fPiv, + 0) — (continuity equation) 



dP 

dv 



P « l^j^j {p{v,t)) 



(Darcy's law) 
(a special type of barotropic fluid) 
(30) 



which describe a special type of barotropic fluid 15]. "A 
special type" means as follows. From the definition of re- 
logarithm ([2]), the last relation can be written as 



Pcx(p^-p-^), 



(31) 
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i.e., a pressure P of this type of barotropic fluid depends 
on the difference between one term with a polytropic index 
K and the other term with —k. 



3.2 Classical Lie symmetry analysis 

In the foUowing we classify the Lie symmetries of the k- 
diffusive equation. 

For the sake of comparison, we first consider the PME 
(|29|) . In this case, the symmetry generators are given by 
(please refer to [16] for more details) 



,(2) ^ ii 



,(3) 



,(4) 



d_ 

dt ' 
d_ 

dv ' 

d ^ d 
v—+2t— , 
ov at 

d , , d 



(32) 
(33) 
(34) 
(35) 



This means that if p(w, t) is a solution of PME, then 

p^^\v,t)^p{e-%,e-^H) , 
p'^^\v,t)^e-'p[v,e^''-^^H) , 



(36) 
(37) 
(38) 
(39) 



are also solutions. 

Generators ([32]) and ([33]) describe translations invariance 
in X- and i-coordinates, respectively and their group in- 
variant solutions (GIS) are merely the constant functions. 
Generators ([34]) and ([35]) describe scaling invariance of in- 
dependent and dependent variables, respectively. In par- 
ticular, the GIS for the symmetry v^^) ]ja,s a kink-like 
shape and in the q ^ 1 limit, it reduces to the Erf(a;) 
function. 

Further, GIS can be obtained by means of linear combi- 
nation of the vectors v^^\i — 1, 2, 3, 4). For instance, from 
the generator v — w'^^ +a tj'^^) obtain the traveling wave 
solutions which arc physically irrelevant because they di- 
verge. Differently, the most general self-similar generator 
is provided by vector v = w^^^ -I- av'^^\ Although the cor- 
responding GIS can not be given in terms of elementary 
functions, explicit solutions can be obtained for special 
values of the constant a. In particular, for a = — 1 we 
derive the Barenblatt scale invariant solution ITTI 




1,0 







^ (c) 






C- 



1 



\Dtr 



where C is the normalization constant. Note that this so- 
lution is nothing but a g-Gaussian with fj, = 2 — q [17 . 
Next, we study the symmetries for the more complicated 
K-diffusive equation. It is worth noting that the PME is 
characterized by the single-diffusive term p^, which leads 
to the GIS ([39]l . while the K-diffusive equation contains 



Fig. 1. Kink-like self-similar solutions of the K-diffusion equa- 
tion p{x, t) = p{z) with z = /t, for k — 0.5 and several 
boundary conditions; (a) p(0) = 0.2 and p'(0) = 0.13; (b) 
p(0) = 0.2 and p'(0) = 0.2; (c) p(0) = 0.2 and p'(0) = 0.6; (d) 
p(0) = 0.2 and p'(0) = 10. 
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two different powers p^^" and p^^'^ . This double nonlin- 
earity in the diffusive term destroys the scahng invariance 
for the dependent variable p{v,t). 
Thus, the only Lie symmetry generators are 



(41) 
(42) 

2t|, (43) 



so that, if p(w, t) is a solution of the K-generalized diffusive 
equation then 



(1) 


d 




" di ' 


,(2) 


d 




dv ' 


,(3) 


d 




dv 



,(2) 
^(3) 



:p{e^'v,e~^H) , 



(44) 
(45) 
(46) 



are also solutions. 

Again, the GIS corresponding to the traveling wave gen- 
erator V = v^^"^ + aw^^-* are physically irrelevant because 
merely constants or divergent. Differently, the only self- 
similar generator is v^^^ and the corresponding GIS has 
a kink-like shape as shown in figure [T] These solutions, 
although not normalizable, can describe the evolution of 
an order parameter belonging to the system. However, as 
compared to the case of PME, there is no scaling gener- 
ator like equation ([55)1 and consequently, the K-Gaussian 
is not an exact solution for the K-diffusive equation. 



4 Asymptotical analysis 
4.1 Numerical simulations 

To further study the long time behavior of the K-diffusion, 
we performed numerical simulations against various initial 
probability distribution functions with different shapes. 
We used a numerical method developed in [18] in or- 
der to study the approach to the equilibrium of a given 
Cauchy problem for the one-dimensional generalized dif- 
fusive equation 



d 



expanding support of p{v,t), Gosse and Toscani [18] de- 
veloped an explicit numerical scheme applicable to one- 
dimensional differential equations of the type given in (|47p . 
We here explain the basic points of this numerical scheme. 
First, we introduce the cumulative distribution function 



r{v,t) 



f p(C,t)dee [0,1] , 

J — oo 



(48) 



associated with the probability distribution p{v, t) and its 
inverse relation 



d , . , , 
— r{v,t) =p(v,t) 



(49) 



By using equation (I47|) and integrating by part we obtain 

-rM^J Piijn^J Hp((.tm 



ov \ ov 



(50) 



Next, since r(v,t) is non-decreasing in the v variable, one 
can define its pseudo-inverse function for a given value 
r e [0, 1] as 

w(f , t) = inf {v e M, such that r{v, i) = f } . (51) 
Then, we have 



r {v{r, t), t) = r 



(52) 



Note that f does not depend on time t. In other words, 
at each time t there exists only one point of v(f, t) whose 
cumulative distribution function r{v,t) equals to f. There- 
fore, differentiating both sides of equation (|52]) w.r.t. t we 
obtain 

d , , s . dr dv dr ^ , „, 

J,rivi,,tlt) = - + --=0. (53) 

In the case of dr/dv ^ we derive the time-evolution 
equation for v{f,t) as 



d , ^ d ^ fdv 



(54) 



>(p(i;,t)), p{v,t^Q)=po{v)>Q , 

(47) 

where po(w) G L^{R) is an initial localized state and «?(p) £ 
C^(M"'"). For computational simplicity, the diffusion coeffi- 
cient D is re-scaled in the independent variables according 
to the symmetry (|43|) . which still holds whatever the func- 
tion <l>{p) is. 

In general, a solutionp(i;, t) of any diffusive equation spreads 
out as time increases. Then, simple conventional methods 
of discretizations in the v-coordinate, which are often em- 
ployed in the numerical simulations, leads to much com- 
putational costs, especially for a long time numerical sim- 
ulation. 

In order to overcome the difficulty due especially to the 



In short, we compute a time-evolution of the pseudo-inverse 
function w(f, t) for a given initial condition f. Since the 
computational domain is now restricted r £ [0, 1] we can 
go around the above expanding support issue for any large 
time. 

For the K-diffusion equation, we set 



•P{p{v,t)) =p{v,t)u^^^{p{v,t)) 



(55) 



The top plot in figure shows a typical time-evolution of 
p(v,t) with an initial state concentrated around u = 0. 
The bottom plot in the same figure [5] shows the time- 
evolution for an initial state with double peaks. Note that 
the double peaks immediately disappear as time evolves 
and become a single peaked shape, which seems to be 
asymptotically approaching to the Ac-Gaussian. 
We have run several simulations with different initial shapes 
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Fig. 2. Top; a typical time evolution for the ^-diffusion equa- 
tion with K — 0.3. Bottom: The same for an initial probability 
distribution function with double peaks. 



p(v,t — 0) to confirm this asymptotic behavior. 

A further indication on this asymptotic behavior can be 

obtained by studying the time evolution of the function 



T]iv,t) = hi 



(56) 



which is the logarithmic ratio between the numerical solu- 
tion p{v,t) of equation (|24p and the trial fitting function 



(57) 



p"(f,i) = exp^^j {a{t)-b{ty) 



This is a K-Gaussian where a{t) and b{t) are the best fitted 
parameters for the numerical solution at each time t. 
In figure [3l we plotted the time evolution of p{v, t) for an 
initial p{v, 0) with a triangle shape. It is clear from this 
picture that ri(v,t) ^ as t ^ oo, i.e. the function (|56p 
gradually decreases to zero as time evolves, which gives a 
strong evidence that the numerical solution is asymptoti- 
cally approaching to the K-Gaussian function. 
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-15000 -10000 -5000 
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Fig. 3. Long time behavior of the numerical solutions with 
an initial triangular probability distribution. The parameter 
K = 0.3 and the number of calculated points for f are 201. The 
quantity r](v,t) of equation (|56p is plotted in each figure. The 
best fitted parameters in (|57p are: a = 1.324, b = 40.54 a± t = 



0.0; a 



-4.862, b = 2.777 x lO"-" at t = 137; a 



-7.640,6 : 



2.871 X 10"* at t = lAx 10^ and a - 



-23.39,6 = 6.205x 10" 



sXt = 6.73 X 10^ , respectively. Inset figures show p{v, t) at each 
time t. 



4.2 Time dependent transformation 

In order to further confirm these asymptotic behaviors 
let us consider a time-dependent transformation for the 
K-generalized FPE (fT9|) . defined by 



w{v,t)^p{m,r{t)) , (58) 
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where 



and 



^W-^, T{t)=\na{t) 
a[t) 



a{t) = (1 + 2ty 



(59) 



(60) 



dt 



w{v,t) = D 



that is the K-difFusive equation with the source term 
w{v,t)/{2t + l). 

In the large time hmit, the source term becomes neghgi- 
ble and equation (|6ip can be weU approximated by equa- 
tion ([M]) . On the other hand, we recall that the solution 
of K-generalized FPE (fT9|) approaches asymptotically to 
the stationary state which is the K-Gaussian ((23|) . There- 
fore, this solution is transformed, by means of equation 
(|58p , into another one that well approximates the solution 
of the K-diffusive equation. In other words, this solution 
asymptotically approaches to a K-Gaussian function. 



5 Summary 



Preliminarily, we observe that the transformation (j58[) maps 8. 
a K-Gaussian into a K-Gaussian. 

By transforming equation according to equation (|58p 9. 
we obtain 

10. 

d , „ 92 /p(v,ty+'^ +p{v,ty^'^\ ^ w{v,t) 

' 2t + l ' 11. 
(61) 12. 



13. 



14. 



15. 



16. 



17. 

We have studied the asymptotic behaviors of the K-diffusive 18- 
equation, which is naturally related with the K-generalized 
statistical physics and the k-FFE (fT9|) . The analysis based 
on the standard Lie symmetry showed that the K-Gaussian 
function is not a scale invariant solution of the K-diffusion 
equation. Notwithstanding, the K-Gaussian plays a rele- 
vant role in the study of the time evolution of localized 
perturbations. In fact, by performing several numerical 
simulations, with different initial shapes, we have found 
strong evidence that these solutions, for large time, are 
well approximated by the K-Gaussian. Arguments based 
on a time-dependent transformation performed on the k- 
generalized FPE also support this result. 
Further confirmations of this behavior could be obtained 
by studying whether a K-Gaussian is invariant under asymp- 
totic symmetries ^lOj of the K-diffusive equation. 
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